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I present here results of a time-dependent analysis of the Dalitz structure of neutral B meson 
decays to 7r + 7r~7r° from a dataset of 346 million BB pairs collected at the T(AS) center of mass 
energy by the BaBar detector at the SLAC PEP-II e + e~ accelerator. No significant CP violation 
effects are observed and 68% confidence interval is derived on the weak angle a to be [75,152] 



I. INTRODUCTION 



when the tagging meson is a B° (B°) is given by 



The time-dependent analysis of the B — > ir + ir~ir° 
Dalitz plot (DP), dominated by the p(770) intermedi- 
ate resonances, extracts simultaneously the strong tran- 
sition amplitudes and the weak interaction phase a = 
zrg[-V td V* b /V ud V: b ] of the Unitarity Triangle [1]. In 
the Standard Model, a non-zero value for a is respon- 
sible for the occurrence of mixing-induced CP violation 
in this decay. p ± ir T is not a CP eigenstate, and four 
flavor-charge configurations (B°(B°) — > /0 ± 7r =F ) must be 
considered. The corresponding isospin analysis [2] is un- 
fruitful with the present statistics since two pentagonal 
amplitude relations with 12 unknowns have to be solved 
(compared to 6 unknowns for the 7r + 7r~ and p + p~ sys- 
tems). 



The differential B° decay width with respect to the 
Mandelstam variables s + , s_ (i.e., the Dalitz plot [3]) 

reads dT(B° -> tt+tt-tt ) = , n i, A l ^ 3 jf ds+ds-, where 
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A 37T {A 37T ) is the Lorentz-invariant amplitude of 



three-body decay _B 



(B° 



V). 



the 
We 



assume in the following that the amplitudes are dom- 
inated by the three resonances p + , p~ and p° and we 



write A 37T = f+A+ + f_A~ + f A° and A 3 



UA + 



f-A~ + JqA , where the f K (with k = {+, — , 0} denoting 
the charge of the p from the decay of the B° meson) are 
functions of s + and s_ that incorporate the kinematic 
and dynamical properties of the B° decay into a (vec- 
tor) p resonance and a (pseudoscalar) pion, and where 
the A K are complex amplitudes that include weak and 
strong transition phases and that are independent of the 
Dalitz variables. 



l^(At)| : 



,-|At|/T B o 

4r R o 



\A 3 „\ 2 + \A 37r \ 2 t (\A 3 „\ 2 - \A 3 ^\ 2 ) cos(Am d Ai) 



± 21m [Astt^L] sin(ATO d Ai) 



. (1) 



where t^o is the mean J5° lifetime and Airid is the B a B° 
oscillation frequency. Here, we have assumed that CP 
violation in b mixing is absent (\q/p\ — 1), Ar^ = 
and CPT is conserved. Inserting the amplitudes A 37T and 
A 3lT one obtains for the terms in Eq. (1) 

L4 3T | 2 ±L4 3 „| 2 = E l^l 2 ^ + 

K e{+,-,o} 

2 E ( Rc U ^ C - Im \f«ffl u ^ m ) - 

K<<re{ + ,-,0} 

Im(A 3w A* 3 „)= E l/«! 2/ « + 

k£{ + ,-,0} 

E ( Re [/„£] j£ + Im [/„£] 1%) , (2) 

K<<re{ + ,-,0} 

The 27 real-valued coefficients defined in Tab. IV that 
multiply the f K fa bilinears are determined by the fit. 
Each of the coefficients is related in a unique way to phys- 
ically more intuitive quantities, such as tree-level and 
penguin-type amplitudes, the angle a, or the quasi-two- 
body CP and dilution parameters [4] (c/. Section IV B). 
We determine the quantities of interest in a subsequent 
least-squares fit to the measured U and I coefficients. 



II. DALITZ MODEL 



With At = t 37r — tt ag defined as the proper time interval 
between the decay of the fully reconstructed B 37T and that 
of the other meson BL the time-dependent decay rate 
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The p resonances are assumed to be the sum of the 
ground state p(770) and the radial excitations p(1450) 
and p(1700), with resonance parameters determined by 
a combined fit to t + — > 77 r 7r + 7r° and e + e~ — > 7r + 7r~ 
data [5]. Since the hadronic environment is different in 
B decays, we cannot rely on this result and therefore de- 
termine the relative p(1450) and p(1700) amplitudes si- 
multaneously with the CP parameters from the fit. Vari- 
ations of the other parameters and possible contributions 
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FIG. 1: Square Dalitz plots for Monte-Carlo generated B° —> 
tv + -k~tt decays. The decays have been simulated without any 
detector effect and the amplitudes A + , A~ and A have all 
been chosen equal to 1 in order to have destructive inter- 
ferences at equal p masses. The main overlap regions be- 
tween the charged and neutral p bands are indicated by the 
hatched areas. Dashed lines in both plots correspond to 
^/s+,_,o = 1-5 GeV/c 2 : the central region of the Dalitz plot 
contains almost no signal event. 

to the B° — > 7r + 7r - 7r decay other than the p's are studied 
as part of the systematic uncertainties (Section IV A). 

Following Ref. [5], the p resonances are parameterized 
in f K by a modified relativistic Breit-Wigner function in- 
troduced by Gounaris and Sakurai (GS) [6]. 

Large variations occurring in small areas of the Dalitz 
plot are very difficult to describe in detail. These re- 
gions are particularly important since this is where the 
interference, and hence our ability to determine the 
strong phases, occurs. We therefore apply the trans- 
formation ds + ds- — ► | det J\ dm' d0', which defines 
the Square Dalitz plot (SDP). The new coordinates are 

to' ee I arccos ( 2 ZiZ m °2* -l), 0' = ±0 O , where 

too is the invariant mass between the charged tracks, 
TO max _ mB0 _ m ^ anc j TO jj" n = 2m w + are the kine- 
matic limits of too and 9q is the p° helicity angle; 9q is 
defined by the angle between the ir + in the p° rest frame 
and the p° flight direction in the B° rest frame. J is 
the Jacobian of the transformation that zooms into the 
kinematic boundaries of the Dalitz plot, shown in Fig.l . 

III. ANALYSIS DESCRIPTION 

The U and / coefficients and the B° — ► 7r + 7r~7r° event 
yield arc determined by a maximum- likelihood fit of the 
signal model to the selected candidate events. Kinematic 
and event shape variables exploiting the characteristic 
properties of the events are used in the fit to discriminate 



signal from background. 



A. Signal and background parametrization 

We reconstruct B° — ► ir + ir~ir° candidates from pairs 
of oppositely-charged tracks, which are required to form 
a good quality vertex, and a 7r° candidate. In order to 
ensure that all events are within the Dalitz plot bound- 
aries, we constrain the three-pion invariant mass to the 
B mass. 

A _B-meson candidate is characterized kinemat- 
ically by the energy-substituted mass toes — 
[(|s + po • Pb) 2 /Eq — p B ]? and energy difference 
AE = Eg - \^/s, where (E B ,p B ) and (E ,p ) are 
the four-vectors of the i?-candidate and the initial 
electron-positron system, respectively. The asterisk 
denotes the Y(4S) frame, and s is the square of the 
invariant mass of the electron-positron system. Wc 
require 5.272 < to es < 5.288 GeV/c 2 . The AE res- 
olution exhibits a dependence on the tt° energy and 
therefore varies across the Dalitz plot. We account 
for this effect by introducing the transformed quantity 
AE' = (2AE - AE+ - AE-)/(AE+ - AE_), with 
AE±(m ) — c± — (ciTcjfmo/mJ 1 ™) 2 , where too is 
strongly correlated with the energy of ir°. We use the val- 
ues c = 0.045 GeV, c_ = -0.140 GeV, c+ = 0.080 GeV, 
m™* = 5.0 GeV, and require -1 < AE' < 1. 

Backgrounds arise primarily from random combina- 
tions in continuum qq events. To enhance discrimination 
between signal and continuum, we use a neural network 
(NN) [7] to combine discriminating topological variables. 

The time difference At is obtained from the measured 
distance between the z positions (along the beam direc- 
tion) of the B®^ and i?° ag decay vertices, and the boost 
/?7 = 0.56 of the e + e~ system: At = Az/ j3~fc. To deter- 
mine the flavor of the -B° ag we use the B flavor tagging 
algorithm of Ref. [8] . This produces six mutually exclu- 
sive tagging categories. 

Events with multiple B candidates passing the full se- 
lection occur in 16% (p ± n T ) and 9% (p°ir°) of the time, 
according to signal MC. If the multiple candidates have 
different 7r° candidates, we choose the B candidate with 
the reconstructed ir° mass closest to the nominal ir° mass; 
in the case that both candidates have the same 7r°, we 
pick the first one. 

The signal efficiency determined from MC simulation 
is 24% for B° -► p ± nT and B° p°n° events, and 11% 
for non-resonant B° — > 7r + 7r~7r events. 

Of the selected signal events, 22% of B° — > p ± n T , 
13% of B° — * p°7T , and 6% of non-resonant events are 
misreconstructed. Misreconstructed events occur when a 
track or neutral cluster from the tagging B is assigned 
to the reconstructed signal candidate. This occurs most 
often for low-momentum tracks and photons and hence 
the misreconstructed events are concentrated in the cor- 
ners of the Dalitz plot. Since these are also the areas 
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where the p resonances overlap strongly, it is important 
to model the misreconstruced events correctly. 

We use MC simulated events to study the background 
from other B decays. More than a hundred channels 
were considered in preliminary studies, of which twenty- 
nine are included in the final likelihood model. For each 
mode, the expected number of selected events is com- 
puted by multiplying the selection efficiency (estimated 
using MC simulated decays) by the world average branch- 
ing fraction (or upper limit), scaled to the dataset lumi- 
nosity (3f fb _1 ). The selected on-resonance data sample 
is assumed to consist of signal, continuum-background 
and B-background components, separated by the flavor 
and tagging category of the tag side B decay. The sig- 
nal likelihood consists of the sum of a correctly recon- 
structed ("truth-matched", TM) component and a mis- 
reconstructed ( "self-cross- feed" , SCF) component. 

B. Dalitz and At distribution 

The Dalitz plot PDFs require as input the Dalitz plot- 
dependent relative selection efficiency, e = e(m',6'), and 
SCF fraction, /scf = fsCF(m' ,9')- Both quantities are 
taken from MC simulation. 

Away from the Dalitz plot corners the efficiency is uni- 
form, while it decreases when approaching the corners, 
where one of the three particles in the final state is close 
to rest so that the acceptance requirements on the par- 
ticle reconstruction become restrictive. Combinatorial 
backgrounds and hence SCF fractions are large in the 
corners of the Dalitz plot due to the presence of soft neu- 
tral clusters and tracks. 

The width of the dominant p(770) resonance is large 
compared to the mass resolution for TM events (about 
8 MeV/c 2 core Gaussian resolution). We therefore neglect 
resolution effects in the TM model. Misreconstructed 
events have a poor mass resolution that strongly varies 
across the Dalitz plot. It is described in the fit by a 
2 x 2-dimcnsional resolution function, convoluted with 
signal Dalitz PDF. 

The At resolution function for signal and B- 
background events is a sum of three Gaussian distribu- 
tions, with parameters determined by a fit to fully recon- 
structed B° decays [8]. 

The Dalitz plot- and At-dcpcndcnt PDFs factorize for 
the charged- B-background modes, but not necessarily for 
the neutral- B background due to B°B° mixing. 

The charged B-background contribution to the likeli- 
hood parametrizes tag- "charge" correlation (represented 
by an effective flavor-tag-versus-Dalitz-coordinate cor- 
relation), and therefore possible direct CP violation in 
these events. 

The Dalitz plot PDFs are obtained from MC simula- 
tion and are described with the use of non-parametric 
functions. The At resolution parameters are determined 
by a fit to fully reconstructed B + decays. 

The neutral-B background is parameterized with 



PDFs that depend on the flavor tag of the event and, 
depending on the final states they can show correla- 
tions between the flavor tag and the Dalitz coordinate. 
The Dalitz plot PDFs are obtained from MC simulation 
and are described with the use of non-parametric func- 
tions. For neutral- B background, the signal At resolution 
model is assumed. 

The Dalitz plot of the continuum events is 
parametrized with an empirical shape, extracted from 
on-resonance events selected in the Toes sidebands and 
corrected for feed-through from B decays. The contin- 
uum At distribution is parameterized as the sum of three 
Gaussian distributions with common mean and three dis- 
tinct widths that scale the At per-event error, all deter- 
mined by the fit. 

IV. RESULTS 

The maximum-likelihood fit results in a B° — > 7r + 7r~7r 
event yield of 1847 ± 69, where the error is statistical 
only. For the U and I coefficients, the results are given 
together with their statistical and systematic errors in 
Table IV. The signal is dominated by B° — > p ± ir T de- 
cays. We observe an excess of p Q Tr a events, which is in 
agreement with our previous upper limit [9] , and the lat- 
est measurement from the Belle collaboration [10]. The 
result for the p(1450) amplitude is in agreement with the 
findings in r and e + e~ decays [5]. For the relative strong 
phase between the p(770) and the p(1450) amplitudes we 
find (171 ± 23)° (statistical error only), which is compat- 
ible with the result from r and e+e~ data. 



A. Systematics studies 

The most important contribution to the systematic un- 
certainty stems from the modeling of the Dalitz plot dy- 
namics for signal. We evaluated this by observing the 
difference between the true values and Monte Carlo fit re- 
sults, in which events are generated based on an alterna- 
tive model. The alternative fit model has, in addition, a 
uniform Dalitz distribution for the non-resonance events 
and possible resonances including / (980), /2(1270), and 
a low mass S'-wave a. The fit does not find significant 
number of any of those decays. However, the inclusion 
of a low mass 7r + 7r _ S'-wave component significantly de- 
grades our ability to identify p 7r° events. . 

We vary the mass and width of the p(770), p(1450), 
and p(1700) within ranges that exceed twice the errors 
found for these parameters in the fits to r and e + e~ 
data [5], and assign the observed differences in the mea- 
sured U and / coefficients as systematic uncertainties. 

To validate the fitting tool, we perform fits on large MC 
samples with the measured proportions of signal, contin- 
uum and B-background events. No significant biases are 
observed in these fits, and the statistical uncertainties on 
the fit parameters are taken as systematic uncertainties 
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'Quasi twobody" Tjf = \A K \ 2 ± \A K \ 



[/ + p l V fit fraction 
Ul p~ii + fit fraction 



0.237 ± 0.053 ± 0.043 
1.33 ±0.11 ±0.04 



Uq Direct CPV (^r 17 ) 
UZ Direct CPV (p~tt + 
U7 Direct CPV (p + tt" 



-0.055 ± 0.098 ±0.13 
-0.30 ±0.15 ±0.03 
0.53 ±0.15 ±0.04 



" Quasi twobody" I K = Im 



A K A K 



I Int. Mixing CPV p'V 
7_ Int. Mixing CPV p~7r+ 
/+ Int. Mixing CPV p+tt" 



-0.028 ± 0.058 ±0.02 
-0.03 ±0.10 ±0.03 
-0.039 ±0.097 ±0.02 



'Interference" L^ Kc(Im) = Re(Im) \A K A°* ± 



,Rc 



ut 

u +~ 

rr+,Re 

U +0 

rr-,Im 

^+0 

rr-,Re 

U +0 

rr+,Im 

u -0 

rr+,Re 

u -l) 

ry-,Im 



0.62 ±0.54 ±0.72 
0.13 ±0.94 ±0.17 
0.38 ±0.55 ±0.28 
2.14 ±0.91 ±0.33 
0.03 ±0.42 ±0.12 
-0.75 ±0.40 ±0.15 
-0.93 ±0.68 ±0.08 
-0.47 ±0.80 ±0.3 
-0.03 ±0.40 ±0.23 
-0.52 ±0.32 ±0.08 
0.24 ±0.61 ±0.2 
-0.42 ±0.73 ±0.28 



" Interference" 


rRe 


= Re 


A K A a * — A" A K * 




r Rc 

1 + - 






-0.1 ±1.9 ±0.3 


T Re 
^+0 






0.2 ±1.1 ±0.4 


rRc 

J-o 






0.92 ±0.91 ±0.4 


" Interference" 


rim 


= Im 


A K A°* + A a A K * 




rim 

J +- 






-1.9 ±1.1 ±0.1 


rlm 






-0.1 ±1.1 ±0.3 


rim 






0.7 ±1.0 ±0.3 



TABLE I: Definitions and results for the 26 U and / observ- 
ables extracted from the fit. We determine the relative values 



of U and / coefficients to Ul 



Another major source of systematic uncertainty is the 
B-background model. The expected event yields from 
the background modes are varied according to the uncer- 
tainties in the measured or estimated branching fractions 
Since £?-background modes may exhibit CP violation, the 
corresponding parameters are varied within appropriate 
uncertainty ranges. 

Continuum Dalitz plot PDF is extrapolated form tues 
sideband, and large samples of off-resonance data with 
loosened requirements on AE and the NN are used to 
compare the distributions of m' and 6' between the ibes 
sideband and the signal region. No significant differences 
are found. We assign as systematic error the effect seen 
when weighting the continuum Dalitz plot PDF by the 
ratio of both data sets. This effect is mostly statistical 
in origin. 

Other systematic effects due to the signal PDFs com- 
prise uncertainties in the PDF parameterization, the 
treatment of misreconstructed events, the tagging per- 
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FIG. 2: Confidence level functions for a. Indicated by the 
dashed horizontal lines are the confidence level (C.L.) values 
corresponding to la and 2<r, respectively. 



formance, and the modeling of the signal contributions 
and are estimated using arious data control samples. 



B. Intepretation of the results 

The U and I coefficients are related to the quasi-two- 
body parameters (Tab. IV B) defined in Ref. [4], explic- 
itly accounting for the presence of interference effects, 
and are thus exact even for a p with finite width. The 
systematic errors are dominated by the uncertainty on 
the CP content of the i?-related backgrounds. One can 
transform the experimentally convenient, namely uncor- 
rected, direct CP-violation parameters C and *4 p7r into 
the physically more intuitive quantities A^~ 



pTT 

and A 



The significance, including systematic uncertainties and 
calculated by using a mininum x 2 method, for the ob- 
servation of non-zero direct CP violation is at the 3.0a 
level. 



C = (C+±C*-)/2 
S = (S+ + S~)/2 
AC = (C+ -C")/2 
AS = (S + - S~)/2 

A - U+ +- U+ - 



A+- 
A~+ 



UT +U 



-+|2_j 



0.154 ±0.090 ±0.037 

0.01 ±0.12 ±0.028 
0.377 ±0.091 ±0.021 
0.06 ±0.13 ±0.029 

-0.142 ±0.041 ±0.015 

0.03 ±0.07 ±0.03 
-0.38 tg;^ ±0.07 



TABLE II: Quasi twobody parameters definition and results, 



and S ± = 



2_r± . 



where C 
and k h = 

diagrams where the p (n) meson is emitted by the W bo- 

A„+C+A„AC 



{q/p)(A+ /A~), so that Aj 



= (q/p)(A-/A+) 
(Ap^~) involves only 



son. A^ and A p ^~ are evaluated as 



l+AC+A p7r C 

Their correlation coefficient is 0.62. 
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The measurement of the resonance interference terms 

allows us to constrain the relative phase <5_| = 

arg (A + *A~ ) between the amplitudes of the decays B° 



p-n+ and B° 



p + ir 



This constraint can be improved 



with the use of strong isospin symmetry. The amplitudes 
A K represent the sum of tree- level (T K ) and penguin- 
type (P K ) amplitudes, which have different CKM fac- 
tors. Here we denote by re the charge conjugate of re, 
where = 0. We define [11] A K = T K e~ ia + P K and 



P K , where the magnitudes of the CKM 



factors have been absorbed in the T K , P K , T K and P K . 
Using strong isospin symmetry and neglecting isospin- 
breaking effects, one can identify P° = —(P + + P~)/2 
and 9 unknowns have to be determined by the fit. 

We find for the solution that is favored by the fit 
<5 + _ = (34 ± 29)°, where the errors include both sta- 
tistical and systematic effects, but only a marginal con- 
straint on <5h is obtained for C.L. < 0.05. 

Finally, following the same procedure, we can also de- 
rive a constraint on a. The resulting C.L. function versus 
a is given in Fig. 2 and includes systematic uncertain- 
ties. Ignoring the mirror solution at a + 180°, we find 
a e (75°, 152°) at 68% C.L. No constraint on a is 
achieved at two sigma and beyond. 



inated by the p resonance. The results are obtained from 
a data sample of 346 million T(AS) — > BB decays. We 
perform a time-dependent Dalitz plot analysis. From the 
measurement of the coefficients of 26 form factor bilin- 
ears we determine the three CP-violating and two CP- 
conserving quasi- two-body parameters, where we find a 
3.0cr evidence of direct CP violation. Taking advantage 
of the interference between the p resonances in the Dalitz 
plot, we derive constraints on the relative strong phase 
between B° decays to p + n~ and p~ir + , and on the an- 
gle a of the Unitarity Triangle. These measurements are 
consistent with the expectation from the CKM fit [12]. 
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